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1. INTRODUCTION 



This is one of a series of reports dealing with modeling of the marine 
atmospheric boundary layer. The basics of the Naval Postgraduate School 
modeling approach are given in Ehirall (1981). That report also describes the 
current status of modeling the pertinent physical processes. The utility of 
the model for tactical use and initial model validation are given in Ehvidson et al 
(1982). This report is in a somewhat different vein. The other reports deal 
with use of the model to describe or predict the evolution of boundary la^er 
properties. Here we use the model "in reverse" to determine an important 
physical parameter, the rate of generation of marine aerosols. 

There are two components to aerosols in the marine boundary layer: 1) con- 
tinental (background) and 2) locally generated sea-spray droplets (Earnhardt 
and Streets , 1970) • The sea-spray droplets are generated primarily by the 
bursting of bubbles at the sea surface (Blanchard and Wbodcock, 1957). The 
bubbles are produced by biological activity, chanical reaction and breaking 
waves (white caps) . Thus, the sea surface is a continuous source of sea-salt 
aerosols in the marine boundary layer. At wind speeds greater than about 3 
meters per second, whitecaps are the primary contributor to the bursting bub- 
bles. These surface produced aerosols can be characterized by a surface flux 
spectrum, F (r), which represents the volune of aerosol per particle radius 
interval produced per square centimeter of ocean surface each second as a func- 
tion of aerosol particle radius. This quantity is a function of wind speed. 

The continuous production of sea-salt aerosols is balanced by several 
removal mechanisms. Che obvious mechanism is the loss of particles as they 
fall back to the surface. This settling under gravity is called "Stokes fall- 
out" and is characterized by the Stokes velocity. The particles are trans- 
ported vertically by turbulence in the marine boundary layer and maintained at 
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a nearly in i form mixing ratio throughout the mixed la^er . The growth of the 
height of the layer constitutes another loss mechanism called "entrainment" 
(Deardorff , 1976) . The final loss mechanisn is "rainout" v^iich occurs vhen 
the particles becane condensation nuclei in the formation of clouds. 

Given the surface flux spectrun and a parameterization of the removal 
mechanians, one could predict evolutions of the aerosol density spectrun 
(Fhirall et al, 1982). Lhfortinately, the surface flux spectrun is not knovn . 
Hawever, v*e can reverse the process — that is, use parameter! zat ions of the 
removal processes and evolutions of the aerosol spectrun — to obtain estimates 
of the surface flux spectrun. This paper describes the calculation of F g (r) 
from data taken by the NPS Ehv ironmenta 1 Fhysics Group during the Cooperative 
Experiment for West Cbast Oceanography and Meteorology (CEWCCM-78) and the 
Joint Air-Sea Interaction Experiment (JASIN) . Details about the experiments, 
equipment, measurements, and analysis are given by Fhirall (1981) • 
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2. AEROSOL SPECTRUM 



The first step to analyzing the evolution of the aerosol spectrun is to 
remove those variations due to changes in the ambient relative hunidity. This 
is dene by transforming each spectrun to a reference relative humidity (RH = 
80% or a saturation ratio S = 0-8) • We make the usual assumption that the 
particle radius at saturation S is given by (Fitzgerald, 1975) 

r g = r*g(S) (1) 

where r is the radius at S = 0.8 and g(S) = 0.81 exp [0.066S/(1.058-S)]. If 
we have a measured aerosol number density spectrum, n' (rg) , then the trans- 
formed spectrum at standard saturation, n(r) , is: 

dn/dr = n(r) = n 1 (r-g(S) ) • g(S) (2) 

Since we prefer to work with the volume spectrum V(r) = 4/3irr n(r) , the cor- 
responding relationship is : 

dv/dr = V(r) = V' (r-g(S))/g 2 (S) (3) 

We now must separate the measured spectrum into its two basic components : 

1) The background continental aerosol of nonlocal origin, V , which is 
present above and within the marine layer; and 

2) The sea-salt aerosol locally generated at the sea surface, v , which 
is well mixed throughout the marine layer (h is the height of the 
mixed la^er) . Thus, we assume: 



V(r) = V c (r) + V s (r) 


Z < h 


(4a) 


V(r) = V c (r) 


Z > h 


(4b) 



The continental aerosol component is represented with a Junge type dis- 
tribution , 

V c (r) = A/r (5) 
where A is the continental coefficient. Since there is very little sea-salt 
volume production for r < 0- 3 urn, we can use the small size range of the 
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spec trim to calculate the continental coefficient (A = rV(r)) . This is docu- 
mented in Fig . 1 where ensemble averages of volume spectra indicated that V(r) 
is essentially independent of wind speed at r = 0*1 pm. At low wind speeds, 
the r = 0.1 urn and r = 0.3 pm ensemble average spectra are described quite 
accurately by Eqn (5) . Ihe claim that A is an accurate index of continental 
influence is nicely validated by comparison with atmospheric fed on activity 
(Iars^n e t al, 1979) from CEWCCM-78 (Fig. 2) • Thus, the locally generated 
sea-salt component can be calculated using: 

V s (r) = V(r) - A/r (6) 

The remainder of this report will deal exclusively with the volume spectra 
transformed to standard humidity, V (r) . We write this as V to simplify the 

o 5 

notation — with the understanding that V refers to the spectrum at RH = 80%, 
with radius, r, an implied variable. 

The evolution of the aerosol spectrum in the well mixed marine boundary is 
described by the following equation (feirall et al, 1982), 

hdV s /dt = F s - (W e + W s ) V s (7) 

( which is presented here without further discussion) 

vhere W g is the entrainment velocity and W g is the Stole s velocity ( Wj , 1979) 

W g = 1.57 x 10 -2 (1 + 1.2(0.8l/g(S) ) 3 )r 2 g 2 (S) (8a) 

W g = dh/dt - W (8b) 

vhere r is in pm, W g is in am/ sec and W is the mean vertical air motion due to 
synoptic scale weather (subsidence) . In Eqn (7) we have left out the cloiri 
formation removal mechanism because it acts on a much longer time scale . Che 
assumption implicit in Eqn (7) is that the convective mixing velocity, W* 
(feimal et al, 1976), is much larger than W g . Since W* is typically 1 m/sec 
vhile W g is on the order of 0. 1 m/sec for the largest particles ( r = 15 pm) 
considered here, the assumption is reasonable. 
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dV/dr (#|4m 2 /cm 3 ) 




Radius (pm) 



Figure 1. Ensemble average total aerosol volume spectra from JASIN. 

The number to the right of the spectrum is the wind speed 
category in m/sec. 
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CEWCOM - 78 




DAY (MAY, 19781 

222 

Figure 2. Time series of atmospheric Rn activity (solid line, Larsen, 
Kaserair and Bressan, 1979) and continental aerosol coefficient, 
A (dashed line) . 
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3. ELUX SPECTRUM ANALYSIS 



The CEWCCM-78 data lave been used to evaluate the surface flux tern of Bqn 
(7). Since the other terms of the expression were measured, one simply calcu- 
lates F g using a rearrangement of the terms. Thus the final expression for 
the flux calculation is, 

F s = Mv s /dt + (W e + W s )V s (9) 

A period was chosen from CEWCCM-78 where all data were available, the 

synoptic conditions were fairly stable, the wind speed was reasonably constant 
and a good mixed layer was present (Thble I) . There was a slowly decreasing 
continental influence indicated by steady 1W winds (A in Table I and Fig. 2). 

Application of Eqn (9) requires a value for the entrairment velocity. If 
the subsidence velocity is known, then W Q can be calculated from dh/dt using 
Eqn (8b) . Since it is difficult to obtain W to 0.1 cm/s accuracy, it was 
decided to estimate W g empirically from the evolution of the mixed-layer tem- 
perature and water vapor density. This was accomplished by simulating the 
analysis period with a dynamic mixed- layer model. After initializing the 
model with the atmospheric parameters at the beginning of the period, the en- 
trainment velocity was adjusted to yield best agreement with the temperature 
and water vapor density at the end of the period. The mean W g for the period 
was 0-35 cm/s. This implies W = -0.3 cm/s for the first 8 hours and W = 0 for 
the last 12 hours. 

The scatter and uncertainty of the aerosol and mixed layer depth measure- 
ments introduces certain trade-offs between averaging times and statistical 
validity. Fbr short averaging times, the random scatter of the individual 
terms is considerable, leading to large variations in the time derivative 
terms. By increasing the averaging time to 4 hours, the F g "signal-to-noise" 
beccmes more reasonable. A time series of the individual terms of Eqn (9) is 
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TABLE I 



. CEWCCM-78 meteorological and aerosol data for the 
time period analyzed for this report. The aerosol 
volume spectra density dv/dr is given at each radius 
r, in ijm. The height of the boundary layer was 
determined from an acoustic sounder. 



DATE 


TIME 


U 

(m/s) 


RH 

(%) 


h 

(m) 


A 


r=. 3 


r=. 8 


r=2 _ r=5 

(jjm 2 /cm 3 ) 


r=10 


r=15 


5/20 


1130 


7.4 


89.4 


420 


7.8 


15.1 


12.5 


9.79 


2.26 


0.801 


0.437 


5/20 


1230 


7.7 


87.6 


430 


5.2 


9. 15 


12.5 


11.9 


2.08 


0.644 


0.327 


5/20 


1330 


8.2 


86.4 


430 


4.5 


8.39 


11.9 


12.2 


2.91 


1.11 


0. 628 


5/20 


1430 


8.1 


82.1 


415 


4.7 


8.54 


10.6 


11.4 


2.24 


0.780 


0.420 


5/20 


1530 


8.3 


80.9 


437 


4.0 


7.96 


12.0 


11.2 


1.26 


0.296 


0.127 


5/20 


1630 


8.2 


81.1 


437 


4.0 


7.85 


11.9 


13.6 


2.47 


0.795 


0.411 


5/20 


1730 


9.0 


82.5 


443 


3.6 


7.26 


12.9 


13.4 


2.55 


0-854 


0.455 


5/20 


1830 


10.0 


83.8 


430 


2.6 


6.54 


14.0 


13.2 


3.30 


1.37 


0.819 


5/20 


1930 


10.3 


84.3 


420 


2.4 


6.07 


13.8 


13.2 


3.71 


1.73 


1.11 


5/20 


2030 


9.8 


86.8 


360 


2.1 


6.54 


15.9 


14.2 


4.02 


1.63 


0.958 


5/20 


2130 


10.1 


87.7 


377 


1.9 


6.07 


14.9 


13.2 


3.34 


1.22 


0.678 


5/20 


2230 


9.5 


87.6 


397 


2.7 


6.82 


13.6 


12.8 


3.52 


1.42 


0.833 


5/20 


2330 


9.2 


87.6 


400 


2.7 


6.15 


12.9 


12.7 


3.46 


1.46 


0.891 


5/21 


30 


9.6 


87.3 


420 


3.0 


6.14 


13. 1 


11.8 


2.58 


0.993 


0.569 


5/21 


130 


9.2 


88.3 


490 


2.0 


5.69 


13.5 


12.9 


3.33 


1.26 


0.710 


5/21 


230 


8.3 


87.4 


520 


1.7 


5.29 


13.3 


13.4 


3.44 


1.39 


0.818 


5/21 


330 


8.2 


86.1 


560 


2.7 


6.11 


11.8 


10.9 


2.18 


0.731 


0.390 


5/21 


430 


8.5 


83.1 


597 


2.2 


5.55 


11.6 


9.71 


1.78 


0.659 


0.369 


5/21 


530 


8.2 


81.3 


620 


2.0 


5.33 


10.9 


9.75 


1.85 


0.732 


0.431 


5/21 


630 


7.9 


79.7 


680 


1.8 


5.01 


11.1 


10. 7 


1.96 


0.678 


0.365 
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shown in Fig. 3 for two different particle sizes. The average surface pro- 
duction flux spectrim is shown in Fig. 4. This flux spectrum applies to the 
average surface conditions for the entire 20 tour period (average wind speed 
was about 9 m/sec) . The error bars represent the uncertainty in the mean 
estimate . 

The relative contributions of entrainment, gravitational fallout and 
surface production are nicely illustrated by defining an equivalent surface 
production vertical velocity, W_ = F /V so that in equilibrium (dV /dt = 0) 

w F = w e + w s (10) 

Ebr the average conditions found in the analysis period, the entrainment and 
Stokes terms were roughly equal at a particle radius of 5 pm (Fig. 5) . 
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15 



Figure 3. Time series of the terms of Eqn (9) for 2 ym radius and 15 ym radius 

particles. The dashed line is the average contribution of the surface 
production term and should be the sum of the other three terms. 
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Figure 4. The ensemble average surface flux spectrum, F (r) , for the 
CEWCOM-78 analysis period (average wind speedf 9 m/sec). 
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Velocity (cm /sec 




Figure 5. Ensemble average contribution of entrainment (W^) , Stokes 
fallout (W ) and surface flux W F = F g /V g assuming a state 
of dynamic equilibrium. 
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4. FLUX SPECTRUM AND WIND SPEED 



Because of the reasonably constant wind speeds during the 20 hour CEWCCM- 

78 period, we were able to improve the statistical certainty for the F calcu- 

s 

laticn by combining all the data from the period, assuming that the flux would 
be a reasonable representation for U =*9 m/sec wind speed. Thus, we now have 
the surface flux spectrun at a single wind speed. In order to estimate the 
flux at other wind speeds, we note that the right hand side of Eqn (10) is 
nearly independent of wind speed for equilibrium conditions. Therefore, the 
flux at one wind speed can be related to the flux at other wind speeds if the 
equilibrium volume spectra are known: 

F s< u l> = F S< U 2> v s (u 1>/ v s (u 2> (11 > 

Vfe have available from JASIN a large set of ensemble averages of aerosol 

volume spectra at different wind speeds (Fig. 6 and Thble II). It is a simple 

matter to apply this data to Eqn (11) using the CEWCCM-78 aerosol flux and 

equilibrium spectrum to generate the surface volume flux spectra as a function 

of wind speed (Fig. 7 and Thble III) . At a wind speed of 6.0 m/sec, the flux 

can be integrated in radius space to yield a total sea-salt dry mass produc- 

2 2 
ticn rate of 4.8 mg/m /day which compares well wath the 4.4 mg/m /day esti- 

mated by Blanchard (1963) and 4.2 ng/m /day estimated by Kritz and lyncher 

(1980) . 
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TABLE II. 



9 O 

Equilibrium sea-salt aerosol spectra in ym /an 
as a function of radius (ym) and wind speed (m/sec) . 



r + 


0.8 


2 


5 


10 


15 


u 


6 


2.8 


3.4 


3 


0.3 


0.05 


9 


10 


10 


5 


1 


0.35 


11 


18 


25 


13 


6.4 


4.2 


13 


20 


30 


20 


15 


15 


15 


24 


33 


22 


22 


22 


18 


28 


35 


28 


28 


28 


* 


9 


9 


11 


2 


0.8 


0.4 
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TABLE III. Surface sea-salt aerosol flux, F (r), (ym /cm /sec) 
as a function of particle radius ( ym) and wind 
speed (m/sec) . 



r + 
u 


0.8 


2 


5 


10 


15 


6 


1.3 


1.1 


2.5 


1.0 


0.33 


9 


4.5 


3.1 


4.2 


3.3 


2.3 


11 


8.2 


7.7 


11 


21 


27 


13 


9.1 


9.2 


17 


49 


48 


15 


11 


10 


19 


72 


140 


18 


17 


11 


24 


92 


180 


★ 


9 


4.1 


3.4 


1.7 


2.3 


2.6 



CEWCCM-78 Analysis Period 
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5. NON-EQUILIBRIUM 



Fbr time periods of a few hours, the aerosol spectrum may not be in a 

state of dynamic equilibrium. If we rewrite Eqn (7) in the following form: 

dV /dt + V /t = F /h 
s' s' p s' 

where the time constant, t , is 

r 

T p = V(W e + w s ) (13) 

then we see the analogy of aerosols and a capacitor charged by an applied 
"voltage", F , through a "resistance", (W + W ) ^ • In this analogy, the 

o 6 S 

"capacitance" is h. 

The response time of the aerosol density is a strong function of particle 
2 

radius because W ~ r . Values of t for S = 0.8 and h = 400 m are given in 
s p 

Table IV. 



TABLE IV. The equilibrium time constant for aerosols 
at S = 0. 8, h = 400 and W g = 0.4 cm/ sec. 

r,ym 0.5 1 5 10 15 

t , hours 28 22 11 3.5 0.5 

fr 



The boundary layer mixing time, r m , which represents the time required for 
changes in aerosol density to be evenly distributed throughout the marine 
layer is : 



t = h/w * 

m * 



(14) 



Ebr the CEWCCM-78 analysis period (h = 400 m and W* = 0.6 m/sec) we find T m = 
0.16 hours. Therefore, short term variations in mixing volume (dh/dt) will 
lead to changes in the aerosol density because the production response time is 
much slower than the mixing time. Thus, for time periods on the order of one 
hour, changes in the aerosol density (hdV s /dt) will be highly correlated with 
the mixing volume term (V dh/dt) . This effect, which is particularly notice- 
able for atialler particles (see Fig. 3 and Thble III), is shown in Fig. 8. 



20 



dV/dr - A/r (#(. im 2 /cm 3 ) 




Radius (pm) 



Figure 6. Ensemble average equilibrium sea salt volume spectra from 
JASIN. The number to the right of the spectrum is the 
wind speed category in m/sec. 
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Volume Flux 
(|jm 2 /cm 2 /sec) 




1 10 
Radius (|jm) 

Figure 7. Ensemble average surface flux spectra deduced from Eqn (11) 
and Table II. The number to the right of the spectrum is 
the wind speed in m/sec. 
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- 1.5 - 1.0 - 0.5 0 0.5 1.0 1.5 



-dh/dt (cm/sec) 



Figure 8. Changes in aerosol spectral density (dV /dt) versus changes 

in mixing volume (-dh/dt). This graph illustrates the dominance 
of these terms in Eqn (9) for short time periods. 
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